CALT-68-2075 



Order o^/Sq Correction to the Charged Lepton Spectrum in 

b ciui decays 



Martin Gremm and Iain Stewart 
California Institute of Technology, Pasadena, CA 91125, USA 



Abstract 

We compute the a^/3o part of the two- loop QCD corrections to the charged 
lepton spectrum in 6 ^ civi decays and find them to be about 50% of the first 
order corrections at all lepton energies, except those close to the end point. 
Including these corrections we extract the central values A = 0.33GeV and 
Ai = — 0.17GeV^ for the HQET matrix elements and use them to determine 
the MS h and c quark masses, and \ Vcb\- 
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I. INTRODUCTION 



In the last few years numerous theoretical and experimental studies have focused on the 
electron spectrum in semileptonic inclusive B Xcii^e decays. The electron spectrum from 
free quark decays receives both perturbative and nonperturbative corrections. Knowledge 
of the shape of the spectrum can provide insights into nonperturbative effects in B meson 
decays, and thereby also give some information on the weak mixing angle \Vcb\- In the 
framework of Heavy Quark Effective Theory (HQET) it is possible to show that the quark 
level decay rate is the first term in a power series expansion in the small parameter Aqcd / 

For infinitely heavy quarks the free quark model is an exact description of heavy meson 
physics. At finite quark masses the first few terms in the heavy quark expansion have 
to be taken into account. Expressions for these nonperturbative corrections to the lepton 
spectrum are known to order {AgcD/^b)^ HhH and the 0{as) perturbative corrections to 
the free quark decay were given in [Q. 

The dominant remaining uncertainties are the two-loop corrections to the quark level 
decay rate and the perturbative corrections to the coefficients of the HQET matrix elements 
in the operator product expansion. Here we examine the former. While a full two-loop 
calculation of the electron spectrum is a rather daunting task, it is possible to calculate 
the piece of the two- loop correction that is proportional to f3o = 11 — 2/3nf with relative 
ease by performing the one-loop QCD corrections with a massive gluon. The a^/?o parts 
of the two-loop correction may then be obtained from a dispersion integral over the gluon 
mass If there are no gluons in the tree level graph, the a^/^o part of the two-loop 
contribution is believed to dominate the full result because /5o is rather large. Several 
examples supporting this belief are listed in 

A recent calculation of the a'^Po correction to the total inclusive rate for B —* Xcdvi 
decays showed that the a^/?o parts of the two- loop correction are approximately half as big 
as the one-loop contribution, resulting in a rather low BLM scale of hblm = O.lSmb. For 
the electron spectrum we find that this part of the second order correction also amounts to 
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about 50% of the order as contribution, at all electron energies except those close to the 
endpoint. Close to the endpoint the corrections are roughly equal in magnitude. 

The HQET matrix elements Ai and A can be extracted from the electron spectrum in 
B Xcii^e decays |T0[. Even though this method of obtaining HQET matrix elements 
was found to be rather insensitive to the first order perturbative corrections, it is useful 
to extract Ai,A including the a^/?o corrections. Then these matrix elements can be used 
to relate the pole quark mass to the MS masses at order a'^Po- Similarly, one can include 
the ttf/3o parts of the two-loop contribution in the theoretical prediction for the total rate, 
which is needed for the determination of |V"cb|. Since the a^jSo corrections are rather large 
the resulting changes in the quark masses and \Vcb\ are not negligible. 

In Sect. II and III we give analytic expressions for the contributions from virtual and 
real gluon radiation. The last phase space integral in the virtual correction and the last two 
integrals in the bremsstrahlung are done numerically. Readers not interested in calculational 
details are advised to skip these sections. In Sect. IV we combine the results from the 
previous two sections to obtain the a'j.Po corrections to the electron spectrum, and discuss 
the imphcations for the extraction of A,Ai, the MS quark masses, and \Vcb\- In Appendix 
we give an interpolating polynomial which reproduces the two-loop correction calculated 
here. 



II. VIRTUAL CORRECTIONS 

The corrections from massive virtual gluons can be calculated in complete analogy to 
the usual one-loop QCD corrections. The ultraviolet divergence in the vertex correction 
cancels when combined with the quark wave function renormalizations. There is no infrared 
divergence since we do the calculation with a massive gluon. The virtual one-loop correction 
to the differential rate can be written as 

dy " ' 48^^^ J ~ 9 )w + 1 - r -y){ai + a^r) - ^rq aa 

+{q\y-i) + y{i-r')-y>3\ (i) 



where fi = fi/ nib is the rescaled gluon mass, and y = 2Ee/mf,, r = mc/mi,, and = q'^/ml 
are the rescaled electron energy, charm mass, and momentum transfer, respectively. The 
limits for the integration over are 



< g2 < 



y{l-y- r^) 



0<y <l-r^. 



(2) 



The functions att,r(g^) and ai(g^), i = 1,2,3 are the contributions from the wave function 
renormalization and the vertex correction respectively. They can be expressed in terms of the 
scalar two- and three-point functions Bq and Co fll]], and the derivative Bq = dBo{a, b, c) / da. 
Explicit expressions for these functions are given in Appendix 0. Using the standard de- 
composition for the vector and tensor loop integrals [0] we obtain 



ai = -2 + 4Coo + 2{Cu + Ci + r^Csa + r'Ca) + 2(1 - + r^){Ci2 + Co + Ci + C2) 
a2 = 2r(Ci + C2), as = -4(Cn + C12 + Ci) - 4:r\Cu + C22 + C2), 

^ 2 - fio(l, A', 1) - 5o(r^ /i^ r^) + (1 - (5o(l, /i^ 1) - fio(0, /i^ 1) 



(r2 - 



5o(r^ /i^ r^) - 5o(0, /i^ r^) + 2(2 + /i^)5^(l, /i^ 1) 



+2(2r^ + /i^)i?^,(r^/i^r^) 



Defining /i = 1 + — and /2 = (/f — 4r^) the coefficient functions take the form 



/2 + /i'(/i - 2)i?o(l, 1, A^) + /i^(/i - 2r^)i?o(r2, r^, /i^^ 



Cn = - + 

/2 



+ (/2 + 2gV)5o(r , 1, r') + 2/i^(/2 + g^ )Co(l, g^ fi\ 1, r^) 



2g^/2 

3r2/i2(/i - 2r2) 



2/2 



-5o(0,l,/i^ 



2 2 ^2^ , 2gV(/2 + 6gV)-/2(/2 + 2gV^ 



(3) 

(4) 



(5) 
(6) 



5o(r^r^/i^) + 



2gVI 



5o(g^l,r^ 



A-(6rU-2)-/2(A + 2)) 
H 2^ ^0(1,1,/^ j 

, 2/iV/2 + /i^(/2 + 6gV) ^ .2 2^2. 2, 
+ Co{l,q ,r ,/i ,l,r ), 



G 



fi 

1 , (l-r2)(2-/i) 
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/2 



2gV2 



5o(0,l,r2) + 



2g^A2(/2 + 6g^)-/2(/2 + 2g2 

2gVI 



-i?o(gM,r^ 



H 2r^2 — ^ ' J J2 ^oll, 1, J 

A2(6r2(/,-2r2)-/2(/i + 2r 



2r2/| 



/i2(2/2 + /i2(/2 + 6a)^,^ 

Cn 1, 



(l,g^r^/i^l,r^), 



Ci2 = H ^TTT Bo{0, 1, r ) H So(0, r , /i J 



2/2 



2g2/2 



/2 



A-(6(A-2r^)-/.) A-(6rU-2)-/2) 
^772 -DoU, -L,/^ jH ^772 -°o 



2/1 

/2(/2 + gVi) - 2AY(3gVi + /2) o . --2 



2/1 



' 2 2 -2\ 



f/l 



i?o(r,l,r^) 



/i'(-/i/2-/i^(3gVi + /2))^ , 

72 ^0(.l; 



/.^ 



"2 Z "2 1 2\ 

q ,r ,/x ,l,r ), 



1 



/i5o(l, 1, fi'') + {2r' - h)Bo{q-\ 1, r') - 2r'B,{r\ r\ fi''' 



Co 



1 



/2 



- 25o(l, 1, A') + (2 - h)B^{q\ 1, r^) + hB,{r\ t\ fi^' 



+/i^(2-/i)Co(l,g^r^/i^l,r^) 



(7) 



(8) 



(9) 



(10) 



The infinite parts of tlie regularized two-point functions can be shown to cancel in eq. 
In the limit fi —>■ the vertex correction diverges logarithmically. This divergence will be 
canceled by corresponding divergences in the bremsstrahlung contributions discussed in the 
next section. 

III. BREMSSTRAHLUNG 



The bremsstrahlung correction is found in the usual manner, by inserting a real gluon 
on the c and b quark lines. The calculation here is complicated by the four-body phase 
space with two massive final states. We follow the standard procedure of decomposing 
the four-body phase space into a two- and a three-body phase space by introducing the 
four- momentum P = Pc + Pg. In the rest frame of the b quark this decomposition reads 



dR4 = dP"^ dR3{mb;pe,Pi?, P)dR2{P;pc,Pg)- (11) 



The 0{as) bremsstrahlung correction to the differential rate is given in terms of dimension- 
less variables (po = P^/rub, = P^/ml) by 
^(1) r rA n'i it/. 12 ^5 

rem, 

dy 

+62(2 - 2P0 - y)y + 63(1 - y - P2)(2P0 y - p-^ - \) 
+64(1 - y - P')y + &5(2P0 + y - 2)(1 - 2P0 - y + P^) 



^rirem^(A) ^ (V) /" ,^2 .po fpO^ _ p2\-5/2 



26i(l-2PO + P^) 



For convenience the above rate has been written in terms of the coefficients 6, 



63 = (PO' - P')ci + (P2 + 2P0)C2 + 3C3 - 3P0(C4 + C5), (15) 

64 = -po{po^ - P^)ci - 3P0p^c2 - SPOca + (po^ + 2p2)c4 + 3P0^C5, (16) 

65 = -po{po^ - P'')ci - 3P0p2c2 - 3P0C3 + 3Po'c4 + (Po' + 2P2)C5, (17) 

which are linear combinations of 



4{vl - vl) 2[{h + 2P0)^ + {ft^ - 2P^f + 2/x^(l + r^)] , 2v+ - fi" 
h ^ h '^^2v_-il^ 



2{vl -vl){2-h) [hfi\3fi^ - 4P0) + 4/x^(2P0 - 1) - IGP" ] 2v+ - fi\ 
h h '^^2v.-fi^' 



+ 



(12) 



61 = (PO' - P2) (P2(C2 - Cl) + Po'ci + C3 - P0(C4 + C5)) , (13) 

62 = (i^°' - P')P'ci + 3P^C2 + (P' + 2P0')C3 - 3P0P='(C4 + C5), (14) 



4(2 + /i^)(/. - 2P0)K - V-) 8[(^ + /x^)(P° - /.) + ^P°]K - V-) .^o^ 
(2t;+-/i2)(2^;_-/i2) , I J 



2(A^ - 4)(2P0 - /i2)(t;^ - v_) A[h\(i^ - P") + 2P°(/x^ + 2P2)](^+ - v_) 

{2v+ - fi^){2v_ - fi^) ' I 

[(P2 + ^2)^/^2 ^ 2(/i2 _ 2P0)(/i - 2P0)) _ ^ 4r2p2;i2] 0^;+ - /i^ 

- h ^^^2^r^^ 

2(2 + /i^)(/2^ - - P^)(2P" - /t^ - /i)(^+ - v_) 
{2v+- ij?){2v_- jx^) 
_ 4[(P2 + r2)(/iP0 - /^^ + /i2i3o) + 4r2p2po](^^ _ 

_ 4(P0-P2)(^2 _^2) 2(2 + /t2)(2P0 - /x^)(/i2 - 2P0 + /t)(^;+ - v^) 
h ^ {2v^ - p?){2v- - fP) 
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(20) 



2[4P0^(2P2 + ^ _ r^)h{2P^ + h- 4P0) + hP\P^ + - 8P^)]{v+ - 

^[2(/i2 - p0)(/,^2 _ 2p0h + 4po2) + /,2^^ _ po + ^2^ _ 2^2 ^^2 ^ ^2po + pO) 

•2t;_ -/i2^' 



+ /i^(l+r^-2P0)]ln(^-^^), (21) 
2(/i4-4)(p2_^r2-^2)(^^_^_) 



_ 2[fi^h^ + (P2 + r2)(2^2 _^ - /^2) + 8r^P^]{v+ - v^) 

2[i,^h - (P2 + ^2) (^2;, ^ 4po) ^ 2PV] , , 2^+ - 

/I ^^^^ 

In the expressions for the q we have put h = — r"^ and 



(p2 + /i2-r2)PO± VPO'-PV(^^ + A^-^^)^- / X 

= • (23) 



The integrals in eq. ([1^) are done numerically between the kinematic limits 

^ ^ <P°<^^^, (24) 



2(1-1/) - - 2 

(/X + r)2 < p2< 1 - y. (25) 

To improve the numerical stability for small /t^ we found it useful to do the P° integral with 
the variable ln(P'' — r^). The remaining limits for this four-body decay are 



0< f^< Jl-y-r, 0<y<l-r. (26) 



IV. THE al(3o CORRECTION 



Combining the corrections from virtual and real gluon radiation, eqs. (|l],|T2D, we obtain 

(1) frA ^r(i) 



ay ay ay ^ 

In the /t ^ limit eq. ( p?!) yields the one-loop correction to the electron spectrum. We have 
checked that our expression reproduces the result in in this limit. The a^Po part of the 
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two-loop correction is related to the one-loop expression calculated with a massive gluon by 
i 

dr(^) ^ _a£2A) 7d^ /dr«(^ _ 1 drW(o) \ 

dy 47r y /t^ y dy 1 + /t^ dy / ■ 

Note that a^^-*, defined in the V-scheme of BLM is related to the more familiar a<j defined 
in the MS scheme by 

a^P = as + l^/3o + ---. (29) 

as is evaluated at unless stated otherwise. In the MS scheme the a^/3o part of the 
two-loop correction reads 

dr(^) ^ 5 a^dr«(0) _ aJ3o 1 df^ / drW(/i) _ 1 drW(0) \ 
dy 3 47r dy 47r 7 /i2 1^ dy I + fi^ dy j ' ^ ' 

The dispersion integral has to be done with some care. We found that using ln(/i2) instead 
of fi^ as the integration variable simplifies the numerical evaluation considerably. 

In Fig. 1 we plot the a^/^o part of the two-loop correction^ and for comparison the one- 
loop correction to the electron spectrum, using r = 0.29, dig = 0.2, nj = 3, and dividing by 
To = G'|'|Kfepm^/1927r^. Except for electron energies close to the endpoint, the a'^Po cor- 
rections are about half as big as the first order corrections. The perturbation series appears 
to be controlled but the higher order corrections clearly are not negligible. Integrating over 
the electron energy we reproduce the result for the correction to the total rate given in Ref. 



In |10] the HQET matrix elements A, Ai were extracted from the lepton spectrum using 



the experimentally accessible observables 



^ J,,,Q,ydE,E,dT/dE, ^ J,,,^^ydE,dT/dE, 

' J,,,^,ydEedT/dE, ' ' A.5Gevdi?,dr/dE/ ^ ^ 



* In Appendix ^ we give a polynomial fit to our results that reproduces the numerical results to 
better than 1% for mass ratios in the range 0.29 < r < 0.37. 
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It is straightforward to calculate the a^Po corrections to these quantities. In the spirit of 
HQET, we use the spin averaged meson masses = 5.314GeV, and rriD = 1.975GeV 
instead of quark masses. Keeping only two-loop corrections that are proportional to /5o, and 
neglecting terms of order a'^Po-^Qcn/'niB we find 

Ri = 1.8059 - 0.035— - 0.082^^ + • • • , (32) 
R2 = 0.6581 - 0.039— - 0.098^^ + • • • , 

where the ellipsis denote the other contributions including nonperturbative corrections dis- 
cussed in []5|,p!0[|. The BLM scales for these quantities are iiblm^Ri) = O.Olm^, and 
fi'BLM{R2) = O.OOTm^, reflecting the fact that the second order corrections are larger than 
the first order. This is a result of the almost complete cancellation of the first order per- 
turbative corrections from the denominators and numerators in Ri^2- In eq.(p^) the BLM 
scales for the numerators and denominators are separately comparable to the BLM scale for 
the total rate Hblm ~ O.lm^. Therefore the very low BLM scales of Ri^2 do not necessarily 
indicate badly behaved perturbative series. 

In order to demonstrate the impact of the a^/^o corrections on the extraction of A, Ai, 
we repeat the analysis of [|l^ neglecting nonperturbative corrections of order {^Qco/^bY 



which may be substantial |^. Because of the higher order nonperturbative corrections, 
large theoretical uncertainties have to be assigned to the extracted values of A, Ai. We 
find that the central values are moved from A = 0.39 ± O.llGeV, Ai = -0.19 ± O.lOGeV^ 
to A = 0.33GeV, Ai = —0.17. The shift in the values of the HQET matrix elements lies 
well within the Icr statistical error of the previously extracted values. However, using the 
values of the HQET matrix elements extracted at a given order in to predict physical 
observables at the same order in a^, guarantees that the renormalon ambiguity in A and 
Ai will cancel []T2|,0] if the expansion is continued to sufficiently high orders in a^. Thus 



including the af/5o parts in the determination of A, Ai allows one to calculate the MS quark 
masses consistently at order a^/Jo- To second order in kqcol'^q and to order a^/^o we have 
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mg{mg) = \ TUMeson - A + H j yl — 1.56 — \ ) , (33) 



where mg is the 6 or c quark pole mass and rriMeson is the corresponding spin averaged 
meson mass. With Oisirrii,) = 0.22, dsirric) = 0.39 this yields mb^nih) = 4.16GeV, mc(mc) = 
0.99GeV for the MS quark masses, albeit with large theoretical uncertainties due to the 
effect of the higher order nonperturbative corrections on the extraction of A, Ai [^. The 
value of mbiiTT'b) is in good agreement with lattice calculations rfibimi)) = 4.17 ± 0.06GeV 
and mb(mb) = 4.0 ± O.OlGeV [Q. The weak mixing angle \Vcb\ can be determined by 
comparing the theoretical prediction for the total rate with experimental measurements. 
Including all corrections discussed in ||10[ we find at order a^/?o 

IK>I ^ 0.043 ( ^-'^-f^'J ij^)'". (34) 
\ 0.105 tb I 



V. CONCLUSIONS 

We have calculated the C(a^/3o) corrections to the electron spectrum in 6 — > ciui de- 
cays which turn out to be rather large, about 50% of the one-loop corrections. These 
corrections can be included in the extraction of the HQET matrix elements A, Ai. We 
obtain A = 0.33GeV and Ai = — O.lTGeV^, both somewhat lower than the values 
extracted at 0{as). Using these values and including O(a^/5o) corrections we obtain 
^bif^b) = 4.16GeV, rnc{mc) = 0.99GeV for the MS quark masses and \Vcb\ = 0.043(i?r(i? — >■ 
Xciiyi) /0.105 X 1.55ps/rB)^/^. These results have large theoretical uncertainties due to the 
effect of nonperturbative corrections of order {AQcn/fnbY on the extraction of A, Ai. 
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APPENDIX A: SCALAR TWO- AND THREE-POINT FUNCTIONS 



Here we list expressions for the scalar functions Bq and Cq [|TT| needed for the calculation 



in Section II. With the cut for logarithms along the negative real axis we have 

So a, 6, c = - - _J - / dx In ^ Al 

e 47rA^e ^ Jo fi^ 

= -- Ht-^) + 2 - 1^(4) + ^+ + ^- (A2) 

e 47rA^e ^ /x^ a;+ X- 

where 



= ^\ . (A3) 



C„(l,j^r^AM,^^) ' 



2 



1 + — g2 — 2ar 
Lz2 f — — — ) - Li2 (— + Li2 ( — — — ) - Li2 





— ^4 / 










( 





L%2 + ] -L%2 \ + L%2 



where 





— Z4 




- 1 


\Z2 


- ^6 




- 1 




- ^8 



+U2 — \-Li2[ + — \-U2\ , (A4) 

KZ'i-Z^) \Z-i-ZsJ \Z3-ZgJ \Z-i~Zg,'' 





- Zr, 




Z2 




- Zj 




zz 








-2- 





- 1 




- ^5 




- 1 


\Z2 


- Zj 


(z. 


- 1 



«- ^2 ' ^0- i^^2_q2_2r2a ' ^^^^ 

zq Zo p,^ ± - 

2:1 = 2:0 + a, Z2 = 7:^ r, Z3 = , 2:4,5 = ^ , (A6) 

[1 — a) a 2r^ 



2^6,7 = TTZ^ , 2^8,9 = • (A7) 



f + l-r^± + 1 - r2)2 - Aq^ + ie 2 - fx^ ± ^(2 - fx^y - 4 

The dilogarithms here are defined as 



Lz2{z) = ~ j^dt ''^^ (A8) 



APPENDIX B: INTERPOLATING FUNCTION 

In this appendix we give an interpolation scheme that makes it easy to reproduce our 
a'^Po corrections to the lepton spectrum for b — > ciu^ decays. The interpolation is done as a 
function of y and r with all other parameters left explicit. 
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^ = To (f<'\y,r) + ^/«(l/,r) + ^fHy,r)] • (Bl) 
The well known tree level result is given by the exact equation 



while the first order function, f^^\ can be obtained from Ref. [^]. The interpolating function 
f^"^^ is found by fitting Chebyshev polynomials, T^, to the energy dependence, and quadratic 
polynomials to the mass ratio r. To improve the accuracy of the fit near the endpoint of the 
lepton spectrum the expansion is given in terms of the variable 



ln(r2) 

Our fit is accurate to better than 1% for the region 



(B3) 



0.09 <?/< 0.99(1 -r^) 0.29 < r < 0.37. (B4) 
The second order function is given by 

f'\y,r) = E f E rA T„_i(y'), (B5) 

n=l \m=0 / 

where Tn{y') = cos(n arccosy'). The 36 coefficients Aij are given in Table 1. Both the 
coefficients as well as Fortran code evaluating dF^^^d?/ and dF^^^dy from the expressions in 
sections II, III, and IV are available from the authors by request (iain@theory.caltech.edu). 
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TABLES 



TABLE I. 


Coefficients An^m for the interpolation 


function f^'^\y,r) 


given in Appendix ^. 


A 


in 


I 


2 


Tt 1 


—4 q 


1 Q 3 


—90 4 


2 


—2.72 


12.4 


— 14 7 


O 




— 19 4 


11 ^1 


A 


9 91 

^ . ^ J. 


—in 1 

±\j. ± 


1 1 q 


o 


1 Q7 


— Q 5Q 


1 1 q 


a 


1 1 1 

± . ± J. 


— 5 08 


LI. WU 


7 




448 


—0 974 


8 


-0.336 


1.44 


-1.62 


9 


-0.319 


1.47 


-1.76 


10 


-0.174 


0.818 


-0.994 


11 


-0.0881 


0.417 


-0.509 


12 


-0.0507 


0.247 


-0.309 
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FIGURES 




FIG. 1. Perturbative QCD corrections to the lepton spectrum. The differential rates are given 
in units of Fq, with r = 0.29, = 0.2, and Uf = 3. 
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